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f^^ Abstract. We study the problem of the computation of the effective diffusion constant of 

w t a Brownian particle diffusing in a random potential which is given by a function V(cf>) of a 

S Gaussian field <j>. A self similar renormalization group analysis is applied to a mathematically 

related problem of the effective permeability of a random porous medium from which the diffusion 
constant of the random potential problem can be extracted. This renormalization group approach 
reproduces practically all known exact results in one and two dimensions. The results are 
confronted with numerical simulations and we find that their accuracy is good up to points 
(^ , well beyond the expected perturbative regime. The results obtained are also tentatively applied 

to interacting particle systems without disorder and we obtain expressions for the self-diffusion 
C^ ■ constant in terms of the excess thermodynamic entropy. This result is of a form that has commonly 

been used to fit the self diffusion constant in molecular dynamics simulations. 
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1. Introduction 

In this paper we will consider the late time diffusion constant associated with a Brownian or 
Langevin particle, in D dimensions, advected by a velocity field which is given by the gradient of a 
random potential. The explicit Langevin equation studied is: 

^ = VV(cf>(X))+ V (t) (1) 

where the local potential V is itself a function of a Gaussian random field 4>. In general (up to an 
overall rescaling of time) the potential V can be written and 

V{<j>) = -(3v{4>) (2) 

where (3 is the inverse temperature and v is the physical potential acting on the tracer particle. In 
this formulation n(t) is a Gaussian white noise of zero mean with correlation function 

(r H (t)r ]j (t , )) = 2S ij S(t-t'). (3) 

Therefore when V is not a linear function of 4> the advecting potential is non-Gaussian. The case 
of diffusion in purely Gaussian potentials has been extensively studied in the literature [2 EJ [3] but 
the non-Gaussian case has received much less attention. The Gaussian case has been studied within 
a variety of approximation schemes and among these schemes the most successful has been the self- 
similar renormalization group method which reproduces exact results in one and two dimensions 
and which in addition is in excellent agreement with numerical simulations in three dimensions 
[U El [3] . The case where V is the square of the gradient of a Gaussian potential arises naturally 
in the case of dipoles diffusing in a random electric field (in the limit where the dipole moment 
cquilibriates very quickly in its local field compared to the time-scales over which the diffusion of 
its centre of mass occurs) [4] . The authors of this paper have examined the case 

v = a^- (4) 

in one dimension [5] , where the diffusion constant can be calculated exactly. In this case it can 
be shown that there is a critical temperature at which the diffusion constant vanishes. Below this 
temperature the diffusion is anomalous (more precisely sub-diffusive) and the exponent associated 
with anomalous diffusion can be computed. The exact results of [5] show that the transport 
properties (the exponent associated with the anomalous diffusion) agree with those obtained by 
a straightforward mapping onto a trap model whose trapping time statistics can be deduced from 
the statistics of the field <\> and the Arrhenius law. There have also been some studies of diffusion in 
(non-Gaussian) potentials generated by the potentials due to a distribution of randomly distributed 
particles which interact with the potential via a fixed deterministic potential [51 [7] . 

A naive application of the self-similar renormalization group to this model at the one-loop level 
is not sensitive to the non-Gaussian statistics of the random potential and fails to predict the 
dynamical phase transition associated with the passage from normal to sub-diffusive transport. In 
this paper we reformulate the self-similar renormalization group approach in such away that the 
results for the Gaussian case are unchanged but we reproduce all known exact results in one and 
two dimensions. In addition we show that the approach works well in other cases by comparison 
with numerical simulations, it predicts the dynamical transition, when there is one, and works 
reasonably well outside the perturbative regime. The basis of our analysis relies on the mapping 



Self Similar Renormalization Group Applied to Diffusion in non- Gaussian Potentials 3 

of diffusion in the random potential to diffusion in a medium of random diffusivity, this is mainly 
done as it simplifies the resulting renormalization group flow. Interestingly as a biproduct our 
analysis recovers an approximate results often used in the computation of the effective permeability 
of random porous media El QUI Q3J EG] 

The underlying Gaussian potential we shall study will be assumed to have a short range correlation 
of the from 

(0(x)0(y)) = A(x - y) (5) 

In a translationally invariant and isotropic system the long-time behavior of the mean-squared 
displacement of the process X described by equation ([TJ) is 

(X 2 ) - 2Dn e t (6) 

where K e is thus the long-time effective diffusion constant of the problem. The Fokker-Planck 
equation describing the evolution of the probability density function (pdf ) for X is: 

^ = V-(VP-PW(£)) (7) 

The problem of diffusion in a medium of random local diffusivity k = exp(y(0)) is described by 
the Fokker-Planck equation 

dP 

— = V-exp(y(0))VP. (8) 

The corresponding stochastic differential equation is 
dX 



ijt = cxp (V («£(X))) W (0(X)) + Vexp (V (<t>(X)))r,(t). (9) 

Again under the assumptions of time-translational invariance and isotropy the late-time behaviour 
of a Brownian tracer particle described by the above Fokker-Planck equation is 

(X 2 ) ~ ZDn^t (10) 

where n p is the associated long-time diffusion constant. The effective diffusion constant «4 is also 

the effective permeability of a random porous medium where fluid flow is described by Darcy's law 

(v) 
and if k is interpreted as a local dielectric constant then n e is the effective dielectric constant of the 

medium (see the review [3]). The effective late-time diffusion constants of the two Fokker-Planck 

equations (0 (gradient flow) and IjHJ) (fluctuating diffusivity) are in fact related via 

(11) 



(exp(V(0)))' 

a result which can be shown in a number of different ways 3, 7, 13 . The effective diffusivity can be 
computed from a statics problem. If one considers the Green's function for the random diffusivity 
problem 

V • exp (V((f>)) VG = -<5(x), (12) 
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the effective diffusivity can be read off from the long-distance behavior of the Green's function or 
equivalently the short wave-length behaviour of its Fourier transform: 

which means that on suitably large length scales Green's function reads 

G- 1 ~ -k^V 2 . (14) 

Note we have dropped the disorder average as we shall assume that k/ is self-averaging. 

2. Renormalization Group Approach 

The basic idea of the self-similar renormalization group [U [2] is to average out the short distance 
components of the random field <p down to some wave-length A and then to write an effective 
diffusion equation, of the same structural form, describing the transport at length scales greater 
than 1/A. 

The Gaussian field is decomposed by defining 

/dk 
— — n exp(ik • x)0(k) (15) 

,k|<A \M) 

where A is an upper ultraviolet scale which is initially infinite. We define a slice of the field <f> at 
inverse length scale A 

f dk 

5</>a(x) = / exp(zk ■ x)<f>(k) (16) 

./A-*A|k|<A (2tt)^ 

The self-similar renormalization group process proceeds by integrating out the slice of the field 5<fr 
whilst assuming that the remaining part of the field 0a-<5A can be treated as a constant at the 
length scale A. The correlation function of the slice of the field <5</>a is given by 

(<^A(x)<ty A (y)> = A 5A (x-y) (17) 

and its Fourier transform is given by 

A 5A (k) = A(k)J(|k|,[A-<5A,A]) (18) 

where I is the indicator function 

I(x, A) = 1 if x e A 

= if x i A. __ (19) 

The field 64>a is thus formally of order y/SA. This means that to order <5A one may write at any 
given point 

Hi = (^a) 



S,= ^)-^ 5A (20) 



S D A D - 1 A (A) 

If we apply the self-similar renormalization group hypothesis to the Green's function in equation 
(|12p we expect that after integrating out the random field down to wave number A that on this 
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inverse length scale the running Green's function obeys a similar renormalized equation of the 
form 

V • exp(y A (0A))VG A = -*(x). (21) 

Here G\ denotes the Green's function averaged over modes of </> of modulus superior to A which 
we denote as: 

G A = (G) A (22) 

Now in this equation we write the field <f>\ = 4>a-&a + 5(f) a. In the running equation for G\ we treat 
</>A-5A as approximately constant to obtain 



exp (va(^a-sa) + y (Va'^a^-a) + (V A (H- SA )) 2 )) 



V • (1 + V A (<f>A-5A)S<t>AWG A = -6(x). (23) 
Under these hypotheses we obtain: 

G A = t 9 - r-, (24) 

exp (y A {4> A - SA ) + S -f (V a (^a-sa) + (V a (0a-8a)) 2 )) 

where the Green's function g A is defined via 

V • (1 + V A (4> A _ SA )Sct> A )Vg A . = -S(x) (25) 
Now averaging over the current momentum slice we find that 

Gasa = 7 : r- (26) 



exp 



(Va^asa) + f (V a (H-8a) + (K(<^a-m)) 2 ) 



Note that the average over the slice of the field S(p A can be carried out in the computation of (gA)- 
Now at length scales A the Green's function <?a should behave as 

The equation determining <?a is of the form 

V • (1 + ip)Vg = -S(x) (28) 

where ip is a Gaussian field with correlation function 

^(x)V(y))= J D(x-y). (29) 

Taking the Fourier transform of equation (|28[) yields 

m = £2 " / J^D k ■ ( k - q)^(q)5(k - q). (30) 

This equation can be iterated diagrammatically and can then be averaged to yield a set of Feynman 
diagrams which can be summed in terms of one-particle irreducible diagrams to write: 

^ - k^W (3l) 
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At order SA (which is simply one loop as the momentum in each loop is SA) we find 

f dq ~ N (k • (k - q)) 2 k 2 f dq ~ N 

^ = j w? d ^ - (k- q ) 2 - d y (4^ (q) (32) 

for small |k|. Note that in principle we have introduced higher order derivatives and interactions 
and so the approach is clearly not exact. However we will see that this approach appears to be 
capturing the essential physics of the problem. The correlation function here is given by 

D(q) = l^A-«A)A(A)/(q, [A - SA, A]) (33) 

and thus we find that 

K (<PA-«a) = 1 J) M ' *- 3 > 

This yields 

Ga-sa ~ t r- (35) 

k 2 exp (Va&asa) + % {V a (4>a-sa) + (1 - %W' 2 a(cPa-sa))) 

we now associate the prefactor of the term k 2 in the denominator above as the effective diffusion 
constant in a region of size 1/A which we denote as exp (Va(4>a-sa)- We now compute the flow of 
the function Va to obtain 

The boundary conditions on Va is V^ = V and the effective diffusion constant is given as 

«W=exp(V (0)) (37) 

i.e. after all the random modes have been integrated out. The renormalization group flow equation 
is non-linear but in the case where V = a<p 2 + b<fi + c the flow does not introduce new interactions 
and the full solution can be computed. However one may formally compute the effective diffusion 
constant via the following observation. If one wants to compute the average 

A = (exp (aU((j>))) (38) 

one may also use a (albeit very simple) renormalization group procedure writing 
A A = (exp (aU(<f>))) A 

= exp(aC/ A (0A)). (39) 

The flow equation for Ua is easy to compute and is given by 

dU A Idn / ,, 2 \ , . 

= oTT \ U A+ aU A) ( 40 ) 



dA 2dA 
Thus if we make the following identification 

V 00 = U 00 ; a = 1 - -^ (41) 

we find that Va = Ua for all A and consequently 



#) = cxp(y o (0)) = exp ^1 - y)j K.(0)) 
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eXP (( 1 -^)^ ) )) F * T " (42) 

Therefore if one has a local diffusivity or permeability k(x) = k(0(x)) i.e. that is an arbitrary 
(positive) function of a Gaussian field, then we find the effective permeability is given by: 

i 

^^ta 1 "^ 1 " 4 * (43) 



This result is a widely used approximation in the field of effective permeabilities and this form is 
sometimes referred to as the Landau-Lifshitz-Matheron conjecture [H [IT] (although is is usually 
stated for Gaussian fields in terms of the local field variance) . This formula is exact in one dimension 
and is also exact in two dimensions if the local permeability is (up to a constant multiplicative factor) 
statistically identical to its inverse [31 HU [15] . If one repeats the argument above for a system 
where 

k = exp(V(cf>) - V{4>')) (44) 

where <f> has the same statistics as <j)' then we find that nf = 1 which agrees with the exact result 
in two dimensions. The result is not exact for the Gaussian case V(4>) — (j> in three dimensions 
but the deviation from the real result in fact only shows up at three loop order |16j . In the 
hydrology community the question whether the Landau-Lifshitz-Matheron conjecture was exact in 
three dimensions animated debate for sometime. 

Now we return to the problem of diffusion advected by the gradient of random potential 
V ((/>), putting together the results of approximate equation ([42]) and the exact relation (fTT|) we 
obtain 



(exp((i-a)y(fl))F3J 

(expTO))) • (45) 

This is the main result of our paper and in what follows we shall analyse the behaviour for 
various choices of the potential V(</>) and confront the predictions with results of numerical 
simulations. 



3. Discussion and some special cases 

In the case of a purely Gaussian potential V(4>) = —/3<f> this gives 

K e = cxpf-— j, (46) 

where we have set the variance of the Gaussian field A(0) = 1. This results is in agreement with 
the renormalization group approaches in references [HE]- It is known to be exact in one and two 
dimensions and it is correct at two-loop order in perturbation theory. However it has been shown to 
break down at three loop order in three dimensions |16] and so the result equation (j45|) is certainly 
not exact. However numerical simulations in three dimensions have shown that the prediction (|46[) 
is remarkably accurate well beyond the perturbative region (where basic perturbation theory should 
work well). 
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A case recently studied by the authors is that where we take 

V{#) = -^, (47) 

i.e. the physical potential v{4>) — </> 2 /2. Unlike the Gaussian case the diffusive behavior will 
depend on the sign of the inverse temperature (3. When f3 is positive the Langevin particle is 
attracted to regions where = 0. In a D dimensional space the regions where (f> — form D — 1 
dimensional subspaces, at low temperatures one thus expects the particle to be confined to these 
regions. However there is no clear mechanism for confining the particle and thus we expect the 
transport to be diffusive at all finite temperatures. When f3 is negative the particle is attracted to 
points where the field (f> is maximal or minimal where V</> = 0. 

In the generic case where attractive regions are zero dimensional they correspond to localised traps. 
In this case the average time to escape from a trap is given by the Arrhenius law 

t~ Tq exp(/3A£) (48) 

where AE is the energy barrier associated with the trap and To a microscopic time scale. Now we will 
assume that effective energy barriers scale like the potential V itself and hence write AE = v— v{4>) 
where V represents an arbitrary energy level at which one is deemed to be not trapped. This gives 
the mean residence time of a trap averaged over traps to be 

(r) ~ to exp(/377) (exp(-/3i/)> (49) 

It is now easy to see that the average number of jumps from trap to trap is given by n ~ 1/(t) and 
thus we find that 

11 

Ke ~ r ~ <exp(-/3,)> (50) 

In terms of the physical potential v our main result equation (|45p reads 



_ < ffl q,(-/3(l-^)^)))^*) ,„. 

(exp(-/M0))) ' [ ' 

Comparing equations (|50|l and (|5ip we see the appearance on the same average (exp (— /3v(<fi))} in 
the denominator; it is the divergence of this term which is thus responsible for the vanishing of the 
diffusion constant. Let us note that, even though this term may diverge at a certain value /3 = j3 c , 
the term in the numerator which is the effective diffusion constant for the effective permeability 



K 



(P) _ 



^(-■Ml-lW)))^ ( 52 ) 



remains finite beyond this value of (3 and thus the random diffusivity problem can have a finite 
diffusion constant while the gradient flow problem exhibits a vanishing diffusion constant. The 
localisation of a dynamical transition, characterised by a vanishing diffusion constant, via numerical 
simulations is notoriously difficult. First the low but finite value of the diffusion constant as one 
approaches the transition means that one must carry out simulations over long time scales to diffuse 
sufficiently to place one in the steady state (time translationally invariant regime) and in order to 
reach the long time regime of the diffusion process. Furthermore, it was shown by the authors [5], 
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that there are finite size effects if one uses a finite number of modes to simulate the random field a 
la Kraichnan [17] , these effects smooth out the dynamical transition in a similar way to which finite 
size effects affect simulations of critical phenomena. It is for this reason that it is sometimes better 
to simulate the random diffusivity problem, corresponding to the gradient flow problem, and then 
deduce the effective diffusion constant for the gradient flow problem via the exact relation equation 
(TTTT) . We will show a numerical example of the effectiveness of this approach later on. 

A specific example where the effective diffusion constant can be explicitly evaluated is 

K0)-i(0-a) 2 , (53) 



here we find 



■ eX P ~ r, U , ol , 2NWi , m ■ ( 54 ) 



[l + ni - m ^) \ D (l + *(l-&))(l+0. 



In the case D = 1 we find 



^9___/ /? 2 « 2 



^=(l-^)'exp^-^j. (55) 

which is an exact result. For D = 2 we find 

„, = (1 + ffl ^-fM_<f?I_L_). < 56 > 

We note that for D = 1 there is a transition for both positive and negative (3. However in higher 
dimensions there is only a transition predicted for (3 negative (J3 = —1), that is to say when local 
maxima and minima of the field <j> behave as traps. In the case where a — the diffusion constant 
vanishes in a power law fashion reminiscent to that predicted by mode coupling type theories [18] . 
In the case where a/0 the dominant behavior in the vanishing of the diffusion constant (or the 
divergence in the characteristic time scale) has the form 

Ke ~ eXP (-^^) (57) 

which has the Vogel-Fulcher-Tammann form often evoked in the analysis of the experimental glass 
transition. Let us note that in one dimension the result for the diffusion constant is a function of f3 2 
and there is a dynamical transition at /3 = 1 and — 1. The transition at /3 = —1 is of course already 
guaranteed due to our formulation of the problem via the diffusivity representation; the transition 
at P = 1 is however predicted directly by the renormalization group analysis. 

The behaviour of K e at positive [3 in two and higher dimensions is very interesting. Recall here the 
particle will be localised on the D — 1 dimensional surfaces cf> = a at low temperatures. For large j3 
the renormalization group prediction is 



for D > 2 and 



^I eXP (~2^ (l + a2)) ,r>,,) 
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for D = 2. We thus see that the dimension plays a crucial role in the low temperature behaviour 
of the diffusion constant, For D = 2 the dominant effect of temperature on the diffusion constant 
is of the Arrhenius form, implying that the crossing of energy barriers is the major mechanism 
contributing to diffusion. However for more than two dimensions there is a simple power law 
behavior of the diffusion constant at low temperature, implying that the effect of energy barriers 
is somehow marginal. This implies that most transport can be achieved without crossing energy 
barriers which are of order 1 and that the particle manages to diffuse while staying close to the 
surface <f> = a. It is worth remarking here that the way in which the Volger-Fulcher-Tammann 
law arises here is identical to the mechanism arising in the glass model of Vilgis [19]. Here the 
local energy barriers are taken to be of the form Ez, where E is the typical energy barrier due to a 
neighbouring atom or molecule in a network type glass, and z is the local coordination number. The 
energy barrier is thus Ez, and if one assumes that z has a Gaussian distribution about an average 
value zq then one easily finds, via the Arrhenius law, the VFT form for the relaxation time. 



4. Numerical simulations 

In this section we test the predictions of the renormalization group analysis against numerical 
simulations of the Langevin equations ^Q and ([9]). In what follows we will generate the random 
fields using the method due to Kraichnan (3] [17] . In the simulations carried out here we found 
that away from any dynamical transition that the results were not significantly changed in going 
between 64 and 128 modes. The results reported here will be in the majority of cases for 128 
modes. The stochastic differential equations for both the gradient flow and random diffusivity 
problem were integrated using second order Runga-Kutta integration schemes developed in (20] [21] 
and reviewed in [3j. In all simulations the effective diffusion constant for a given realisation of 
the disorder was obtained by fitting the mean squared displacement averaged over 10000 particles 
at late times. The time of the simulation was chosen so that particles had typically diffused ten 
or so correlation lengths of the field. The diffusion constant is determined by a fit of the average 
mean square displacement over the last half of the time of the simulation (to ensure that the mean 
squared displacement is well within the linear regime). In three and higher dimensions the late 
time average mean square displacement is fitted with a simple linear fit of the form At + B and 
in two dimensions a logarithmic correction is used i.e. At + Bln(t). Finally, the average over the 
disorder induced by the random field was made over 500 realisations of the field. In all simulations 
the Gaussian field was taken to have correlation function 

A(x) =cxp(-ix 2 ) (60) 



and we concentrates on quadratic forms for the potential of the form v — ((f) — a) 2 /2. Firstly we 
carried out numerical simulations of the diffusivity problem described by the Langevin equation 
(0). Shown in figure ([T]) is the numerically measured value for the diffusion constant in two, three 
and four dimensions compared with that given by the renormalization group prediction for the case 
a = 0. We see that the agreement is excellent up to very large values of |/3| showing that the RG 
approach works well outside the expected perturbabitive regime. The RG prediction appears to 
improve as the dimensions of the space is increased. In figure ^ we show the corresponding curves 
for the case a = 0.5, again we see excellent agreement. 
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Figure 1. Analytical calculation of equation (152 I t for k)P for a = in 2(solid line), 3 (dashed 
line) and 4 (dotted line) dimensions, confronted with direct simulation of the random diffusivity 
problem. The simulation results are shown as circles ( 2D), squares (3D) and stars for (4D). We 
have excellent agreement except at negative (3 except for /3 < — 1 in 4D. This is to be expected 
as locally re(x) can become very big which is obviously a problem for simulation. 




Figure 2. Analytical calculation of equation (152 I t for K e p for a = 0.5 in 2 (solid line), 3 (dashed 
line) and 4 (dotted line) dimension, confronted with direct simulation of the random diffusivity 
problem. The circles stand for 2D, squares for 3D and stars for 4D. We have an excellent 

(p) 
agreement except for p < —1 in 4D. In this case, k)T' become indeed very large which force us 

to use smaller time step and thus the simulations are harder. 
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Figure 3. Renormalization group prediction for re e for the case a = in 2 (solid line), 3 
(dashed line) and 4 (dotted line) dimensions, confronted with numerical results deduced via the 
simulation of the random diffusivity problem and using equation JTTJ to estimate K e (2D circles, 
3D squares and AD stars). We see that the results are excellent and improve on increasing the 
spatial dimension. 



The predictions of the RG analysis can also be directly compared with a simulation of the stochastic 
equation {1} and the results are found to be in excellent agreement for small values of \/3\. However 
near the dynamical transition finite size effects play a role. Direct simulation of the gradient flow 
case also requires much longer simulation times to estimate the asymptotic diffusion constant as 
finite time corrections seem to be more importamt. It is clearly much better to simulate the 
stochastic equation ([9]) and then determine the effective diffusion constant for the gradient flow 
problem using the relation equation (jTTJ) . The results using this method for the case a = are 
shown in figure ((3]). We see that the results for 3 and 4 dimensions are in excellent agreement with 
our analytical calculations and that for 2 dimensions the only discrepancy is at large positive values 
of/3. 



5. Conclusion and Discussion 



We have seen that the renormalization group scheme developed in [T] and [5] can be refined to 
take into account non-Gaussian potentials. The new scheme retains the merit of reproducing 
exactly known results in lower dimensions and has good agreement with numerical simulations in 
cases where exact results are not known. The main breakthrough is that this scheme is capable 
of predicting dynamical phase transitions where the self-diffusion constant vanishes, a transition 
analogous to the glass transition. A key point in the analysis was to apply the approximate 
renormalization group not to the problem of diffusion in a random gradient field but to a 
mathematically related problem of diffusion in a random diffusivity field. This renormalization 
group calculation applied to this problem produces a form of the celebrated Landau-Lifshitz- 
Matheron conjecture from the field of random porous media as given in equation (|4"3"|) . 
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Of course the problem we are considering is one with quenched disorder, in glass formers the disorder 
is thought to be somehow self induced. Let us consider for a moment the problem of N Brownian 
particles (of bare diffusion constant 1) interacting via a pairwise potential u. The Langevin equation 
can be written as 

^| = „-/3V„(X) (61) 

Here the potential is given by 

i/(X) = ^u(X 4 -X i ) (62) 

i<j 

i.e. the energy due to the pairwise interaction between particles. The corresponding permeability 
problem thus has 

k(x) = exp(-/3i/(x)). (63) 

Clearly the system is not disordered but we shall treat it as it were and apply the formula (|43[) to 
estimate the self diffusion constant. Firstly we have 

(«(x)) = ^Zn (64) 

where V is the volume of the system, N the number of particles (assumed indiscernable) and Zjy is 
the canonical partition function for the system. Here we have simply replaced the disorder average 
by the spatial average (this is in fact the correct average to make if one looks at the derivation 
of equation (fTTj) [31 [7J [T3] it is replaced by a disorder average by appealing to ergodicity) . If one 
introduces the free energy per particle /(/3) we find that 

(«(x)) = exp (-N(pf(f3) - In(p) + 1)) , (65) 

where p is the particle density p = N/V. Similarly we denote the dimension of the space of the 
diffusivity problem by D = Nd where d is the physical space dimension and find 



^(x) 1 -*)! 77 *! 



exp 



v( //! Mi- 1 




(66) 



We see that the two quantities above have logarithms which are extensive in N but we expect the 
self-diffusion constant to be intensive, remarkably the ratio of the two above quantities is intensive 
and we find to leading order in N that 

K , = oxp (M! m + iM^i). (67) 

Now we use the trivial thermodynamic identity that the entropy per particle is given by s = 
—df/dT = f3 2 df/d(3 and also the fact that at (3 = we have s(0) = — ln(p) + 1, to obtain 

. .,(«(/?) -*(o)r 

exp 



= exp (Js e *(/3)J , (68) 

where s ex is simply the excess entropy per particle with respect to the perfect gas. Thus the 
approximative line of mathematical reasoning we have followed has lead to a quite interesting 
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relation between a dynamic quantity (the diffusion constant) and a thermodynamic quantity (the 
excess entropy per particle). The celebrated Adam-Gibbs relation for glasses relates the relaxation 
time to the configurational entropy [22j , here we have a relation between the self diffusion constant 
and the full entropy and the relationship is quite different from the Adam-Gibbs form. However for 
some time in the chemical physics literature [23, 24, 25, 26 it has been observed that the diffusion 
constant in molecular dynamics simulations, when written in dimensionless form, can often be well 
fitted (in the denser phase) by the expression 

K = Cexp(As ex ). (69) 

Now our computation is for Langevin systems so one would have to course grain a real molecular 
dynamics to a Brownian level to determine kq (which here we have set to 1) for the effective Langevin 
dynamics; we thus cannot reasonably expect to predict the value C without further study. However 
our naive prediction would be A = 2/d. Numerical simulations 23J have revealed that A varies quite 
weakly depending on the species, and it is reported that A rs 0.65 for hard spheres and A « 0.8 for 
Lennard- Jones fluids, both in three dimensions. In another study it was proposed that A = 1 [25] 
and the subject is still debated and studied (see [55]) for a recent review. Here we predict A = 2/3 
which is in intriguing agreement for the hard sphere result ! The only other analytical derivation 
that we are aware of for relations of the type of equation (|68p is via mode coupling theory (and 
thus quite different to that given here) [27] where the effect of mixtures was also included. It will 
be interesting to see if this last application of our method could be refined to treat mixtures and 
also it should be confronted with numerical simulations of Langevin systems [2 8) . 
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